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Abstract. Completing Aronov et al.'s study on zero-discrepancy matrices for digital 
halftoning, we determine all (m, n, k, I) for which it is possible to put mn consecutive 
integers on an m x n board (with wrap-around) so that each k x I region has the 
same sum. For one of the cases where this is impossible, we give a heuristic method 
^ 1 to find a matrix with small discrepancy. 



A semimagic square is a square matrix whose entries are consecutive integers and 
which has equal row and column sums. One way to generalize this millennia-old concept 
is to specify the sums on regions other than rows and columns. Ingenious constructions 
of squares satisfying various sum constraints have been described by many professional 
and amateur mathematicians. While most of them are interested in adding more and 
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more constraints to make their squares impressive, one can generally consider sum 
conditions on any set of regions. 

Aronov et al. [1] took up this problem for square regions: is there an n x n matrix 
with entries 0, . . . , n 2 — 1 such that every k x k region has the same sum? It is 
\ amusing to note that this variant of the classical problem is motivated by an engineering 

question of finding good dither matrices for digital halftoning, a method to approximate 
a continuous-tone image by a binary image for printing (see their paper for details). 
They showed [1, Theorem 1], using what they call constant-gap matrices, that the 
answer is yes if k and n are even or if n is an integer power of k, and no if k and n are 
relatively prime or if k is odd and n is even. We will solve this problem completely by 
determining all (n, k) for which such matrices exist (Section 1). Our construction of the 
matrices is much simpler even for the cases that have already been settled positively. 
We also give counterexamples to Asano et al.'s conjecture on the smallest possible 
discrepancy when n is odd and k = 2 (Section 2). 

Definitions. For a positive integer N, we write [N] = {0, 1, . . . , N—l}. The remainder 
when an integer x is divided by N belongs to [N] and is denoted by x mod N. 

We consider the slightly generalized setting where the matrices and regions are rect- 
angles instead of squares. Let m and n be positive integers. For an m x n matrix D 
and index G [m] x [n], we denote the (i,j)th entry of D by D{i,j). Any set 

R C [m] x [n] of indices is called a region. The sum of the numbers on R is denoted 
by D(R) = J2uj\ eR D(i,j). The discrepancy of D with respect to a set M of regions 
is the difference between the maximum and minimum D(R) as R varies in When 
it is zero, D is said to be ^-uniform. 

The translate of R by (a, b) G Z 2 is denoted by 

(1) R + (a, b) = \ ((i + a) mod m, (j + b) mod n) : {i, j) G R } C [m] x [n]. 
The set of all translates of R is denoted by R = { R + (a, b) : (a, b) G Z 2 }. 
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Figure 1 . This 5x5 table D has discrepancy 8 with respect to [2] x [2] , 
because 44 < D(R) < 52 for every 2x2 region R. 

By an m x n table we mean an m x n matrix in which each element of [run] appears 
exactly once. We are interested in tables with small (or zero) discrepancy with respect 
to [k] x [/], the set of all k-by-l rectangles (Figure 1). 

I. Zero discrepancy 

The greatest common divisor of positive integers x and y is denoted by gcd(x,y). 
The goal of this section is to show the following: 

Theorem 1. Let m, n, k, I be positive integers with k < m and I < n. Let kl = 
gcd(/c,m) and I' = gcd(Z,n). Then there exists a [k] x [l]-uniform m x n table if and 
only if k! and I' are greater than 1 and k'l'imn — 1) is even. 

This is an immediate consequence of the following Lemmas 2 and 3. 

Lemma 2. A [k] x [l]-uniform m x n matrix is [gcd(/c,m)] x [gcd(Z, n)]-uniform. 

Proof. Let D be a [k] x [/]-uniform m x n matrix. We will show that D is [k'\ x [/]- 
uniform, where k' = gcd(/c,m). We get the conclusion of the lemma by repeating the 
same argument with rows and columns switched. 

For each (i, j) e [m] x [n], the regions [k'] x [I] + (i, j) and [k'] x [/] + (i + k, j) have 
the same sum on D, because each of them combined with [k — k'\ x [I] + (i + k',j) 
makes a k x / rectangle. Thus for each (i, j) e [m] x [n], the rectangles 

(2) [k'\ x [I] + (i + qk,j), q G [m/k'\, 

all have the same sum on D. Since k! = gcd(k,m), these m/k' rectangles cover the 
strip [m] x [I] + (0, j) without overlap. Hence, 

(3) %-D([V\x[(\ + (iJ))= D([k'}x[l] + (i + qk,j)) 

qe[m/k>] 

= D([m] x [I] + (0, j)) = D (M x W + ^3))- 

r£[m] 

Since the rightmost side is a constant independent of by [k] x [I] -uniformity, so 
is the leftmost side. Thus D is [k'\ x [/]-uniform. □ 

Lemma 3. Let m and n be positive integers, and let k < m and I < n be their positive 
divisors, respectively. Then there exists a [k] x [I] -uniform m x n table if and only if k 
and I are greater than 1 and kl{mn — 1) is even. 
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(k,l,n) = (2,1,7) 
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(M,n) = (3, 1,7) (M,n) 
Figure 2. Examples of matrices of Lemma 4. 
One direction is a simple generalization of [1, Theorem 1 (b, c)]: 



(3,2,8) 



Proof of the "only if" part of Lemma 3. Let D be a [&] x [I] -uniform m x n table. It 
is easy to see that D(R) = kl(mn — 1)/2 for each R e [k] x [I]. Since D(R) must be an 
integer, the second claim follows. For the first claim, assume k — 1 for contradiction 
(the case I = 1 is similar). Then D([l] x [/]) = £>([1] x [Z] + (0, 1)) and hence D(0, 0) = 
D(0,l), contradicting the assumption that D is a table. □ 

For the converse, we use the building blocks provided by the following lemma: 

Lemma 4. Let k > 1 and I > be integers and let n be a positive multiple of I. If 
kl(n — 1) is even, then there exists a [k] x [l]-uniform k x n matrix in which each row 
is a permutation of [n] . 



Proof. A [k] x [/]-uniform k x n matrix and a [k'\ x [/]-uniform k' x n matrix stacked 
vertically make a [k + k'\ x [/]-uniform (k + k') x n matrix. Also, a [k] x [/]-uniform 
matrix is [k] x [Z'j-uniform for any multiple I' of I. Therefore, it suffices to construct 
the desired matrix P for the cases (k,l) = (2,1), (3,1) and (3,2) (Figure 2). If 
(k, I) — (2, 1), let 

(4) P(0J)=3, P(l,j) = n-l-j. 
If (k, I) = (3, 1), then n is odd by the assumption; let 

(5) P(0,j)=j, P(l,j)= + modn, P(2,j) = (-2j-l)modn. 



If (k,l) = (3,2), let 

(6) P(0,j) = P(1,j) 



+ -(jmod2), 



P(2,j) = n-l-j. 



□ 



It is easy to verify that P is [k] x [/]-uniform in each case. 

Proof of the "if" part of Lemma 3. We may assume without loss of generality that 
l(mn — l) is even. In this case, both kl(n — 1) and l(m/k — 1) are even, so by Lemma 4, 
there are a [k] x [Z]-uniform k x n matrix P whose rows are permutations of [n], and 
a [/] x [l]-uniform I x (m/k) matrix Q whose rows are permutations of [m/k]. Define 
an m x I matrix T by 

(7) T(a,j) = Q(j, [a/k\)k + (a mod k). 



Then T is [k] x [/]-uniform and its columns are permutations of [m]. Define an m x n 
matrix D by 

(8) D(a, b) = P(a mod k, b)m + T(a, b mod I) 
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Figure 3. Construction of D for (k,l,m,n) = (3,2,9,8). 



(Figure 3). Since P and T are [k] x [/]-uniform, so is D. To see that D is a table, 
suppose that D(a,b) = D(a',b'). By (7) and (8) we see that 

{P(a mod k, b) = P{a mod k, b'), 
Q(b mod I, [a/k\) = Q(b' mod I, [a'/k\), 
a mod k = a' mod k. 

Since P's rows are permutations, the first and the third equation imply that b = b' . 
Since Q's rows are permutations, this and the second equation imply that a = a'. □ 

In the above, we constructed the uniform table as a linear combination of two uniform 
matrices with smaller entries. This idea is due to Euler [3] who gave a construction of 
a semimagic square (that is, a ([1] x [n] U [n] x [l])-uniform n x n table) from a pair 
of special ([1] x [n] U [n] x [l])-uniform matrices called Latin squares. 



2. Finding low-discrepancy tables by ranking 

In this section, we confine ourselves, as Asano et al. [2] did, to the case where 
k = I = 2 and m = n. Theorem 1 states that in this case a uniform table exists if 
and only if n is even. For odd n's, they construct a table with discrepancy 2n, and 
conjecture that it is the smallest possible. This is refuted by our Figures 1 and 4. 
Figure 1 was discovered by an exhaustive search. We describe briefly how Figure 4 
was obtained. 

Define /: [0, l] 2 ^ R by f(x,y) = g(x) + g(y), where 



(10) 



9{x) 



I -{Ax -If if x < 1/2, 
l + (4x-3) 2 if x > 1/2 

(Figure 5). Let a, (3 e [0, 1] and define s: [n] 2 ->■ [0, l] 2 by 

'i + a j + (5 s 
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FIGURE 4. A 31 x 31 table whose discrepancy with respect to [2] x [2] is 27. 



Let H be the n x n table whose (i, j)th entry is the rank of f(s(i,j)) (with some 
tie-breaking rule): 



(12) H(i,j) = \{(i',f) E [nf : f(s(i',j')) < f(s(i,j)) or 

(/ (s(if, j')) =f(s(i,j)) and m' + f <ni+j)}\. 



Figure 4 was obtained by this method with n — 31 and (a, (3) = (0.286, 0). 

To see intuitively why D has small discrepancy, note that a 2 x 2 region in D 
corresponds to the region in H (or its translate) shown in Figure 6. These four cells are 
mapped by s to two nearby points (x±e, y) and another two points (x + 1/2, y+l/2±e) 
(the coordinates are modulo 1). Since f(x,y) = —f(x + 1/2, y + 1/2), the sum of the 
values of / at these four points is almost zero. Thus, assuming that taking the ranks 
does not distort the distribution of values too much, we can expect that D has low 
discrepancy. We add the displacement (a, /3) in (11) in order to reduce the chance of 
ties in the ranking which seem to work adversely. 

As Aronov et al. [1] point out, our problem is analogous to a common situation in 
discrete geometry where we try to arrange discrete objects so that they look close to 
some "balanced" continuous distribution. The constraint peculiar to our problem is 



Finally, define the desired matrix D by 
(13) D((i + j) modn, 
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Figure 5. The functions g and /. 



FIGURE 6. A region in H corresponding to a 2 x 2 square in D (for n — 9). 



that we have to use each number in [mn] exactly once. The ranking technique used 
here may be applicable to other problems with this constraint. However, analyzing 
its performance seems to be hard: although our computer experiment for several n's 
suggests that the above method achieves sublinear 2x2 discrepancy we have no proof 
yet. 
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